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Abstract 

The influence of a magnetic field on the broadening of spectral lines and transition arrays in com- 
plex spectra is investigated. The anomalous absorption or emission Zeeman pattern is a superposition 
of many profiles with different relative strengths, shifts, widths, asymmetries and sharpnesses. The 
a and it profiles can be described statistically, using the moments of the Zeeman components. We 
present two statistical modellings: the first one provides a diagnostic of the magnetic field and the 
second one can be used to include the effect of a magnetic field on simulated atomic spectra in an 
approximate way. 

1 Introduction 

The existence of intense magnetic fields (1 - 1000 MG) for white dwarfs was confirmed by the observation 
of the splitting of spectral lines in the visible and UV range. Pulsars and neutron stars, discovered 
through their spectrum in the range of radio-frequencies and X-rays, have an even more intense magnetic 
field. Magnetic fields of a few kG also affect the radiative accelerations in chemically peculiar (CP) stars 
[Alecian & Stift (2002)]. The emission or absorption lines of an atom are modified by the magnetic fields. 
Thanks to the Zeeman effect [Zeeman (1897)], such fields can be detected at large distances, through 
the measured radiation. The methods differ according to the nature of the stars studied, the magnitude 
and the geometry of the magnetic fields and the spectral resolution. Moreover, the variations of the 
magnetic field of stars during their rotation bring some information about their global geometry In 
spectro-polarimetry (see for instance [Donati & Cameron (1997)]), the circular polarization with respect 
to the wavelength is recorded, which enables one to separate the components of the field, parallel or 
perpendicular to the line of sight. For all these applications, a sophisticated theoretical modelling is 
required; however, due to the huge number of electric-dipolar lines in hot plasmas, the exact calculation 
becomes prohibitive, and even useless due to the coalescence of the lines induced by the other physical 
broadenings. In the present work, we propose a statistical model for the simulation of the impact of strong 
magnetic fields on atomic spectra containing many overlapping lines and in which Zeeman patterns are 
unresolved. 

2 Statistical modelling of the intensity and calculation of the 
moments 

The intensity, detected with an angle of observation 9 with respect to the direction of the magnetic field, 
is given by [Godbert-Mouret et al. (2009)]: 

Ie(E)=( 1 -±^)(I + r(E) + IME)) + ^I (E). (1) 
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In the following, an energy level will be denoted by 7J. The quantity J stands for total angular 
momentum and 7, which can be interpreted as the "signature" of the level, represents all the quantum 
numbers, whose coupling leads to J. Each electric-dipole El line 7J — ¥ 7' J' splits into three components 
associated to the selection rules M'=M + q, where the polarization q is equal to for ir components and 
to ±1 for a± components. In the present work, Hanle effect is not taken into account. For a photon 
energy E, the intensity of the q component of the line 7J — » 7' J' (whith energy and strength i? 7 j j7 'j' 
and S 1 j n 'j') reads: 



I q (E) = S M ,w, q *(E - E lM ,j, - fi B B (gyj.M' - g 7 jAf )), (2) 

M,M- 



f J 1 J' V 

where Sm,M',i = c m,M',i x S iM'J' and Cm,m\<i = 3 ( _ M _ q M > J ■ * is tne line profile, 

which takes into account the other broadening mechanisms (Stark, Doppler, etc.). The contribution of 
the magnetic field to the energy of the state 7JM is equal to /isBg^jM, where B is the magnitude of 
the magnetic field, /is is the Bohr magneton, and g 1 j is the Lande factors of level 7 J. The distribution 
Iq(E) can be characterized, in the referential of E^j^iji, by the centered moments of the distribution 
of the energies weighted by the strengths [Mathys & Stenflo (1987)]. The n th — order centered moment 
reads: 

M n}c = E (gyJ'M' - 9-yjM - M[ q] ) " . (3) 

M,M' 

In this expression, is the center-of-gravity of the strength- weighted line energies (relative to 

E^j^iji and in units of fisB): 

M 9l = f (2(<? 7 j + 9 7 'j') +(g y j-g YJ/ )(J-J')(J + J , + l)) = qg e , (4) 

and c i s ^ ne variance (in units of (jibB) 2 ). It is useful to introduce the reduced centered moments 
[Kendall & Stuart (1969)] defined by: 



a lQ] = 



(5) 



The use of a$ instead of allows one to avoid numerical problems due to the occurence of large 
numbers. The first values are — 1, = and = 1. The distribution I q (E) is therefore fully 
characterized by the values of anc ^ °f the high-order moments a$ with n > 2. The first four 

moments are often sufficient to capture the global shape of the distribution [Pain et al. (2009)] and 0$ 
and aj^ (see table [T]), named skewness and kurtosis, quantify respectively the asymmetry and sharpness 
of the distribution (a^ = 3 for a Gaussian). The moments can be easily derived using Racah algebra 
[Landi DegPInnocenti (1985), Mathys & Stenflo (1987), Bauche & Oreg (1988)]. 



3 Taylor series expansion and diagnostic of the magnetic field 

Assuming, for simplicity, that ^ is a Gaussian with a variance v (which is exact only for Doppler broad- 
ening), the quantity I q (E) of equation @ can be expressed [Mathys & Stenflo (1987)] as a second-order 
Taylor series (Rodrigues' formula): 
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Table 1: Values of Q3 and 0:4 of the <r+ anf 7r components of El lines (0:3=0 for the ir component since 
it is symmetric). 
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The Taylor-Series method is valid for \ibB < y/v, but breaks down if [IbB becomes much larger than 
y/v. However, assuming the knowledge of the variance v of the other broadening mechanisms, expression 
([5]) enables one to estimate the magnitude of the magnetic field from the measurement of the full width 
at half maximum (FWHM) of the line: 



B = — 



v(l - 2e- 52 /2) 



Hb\ C(#)(l-2e- 52 /2(i_,52))' 



where (5=FWHM/(2^), C{6) = A{6) [M [ ^ 
B{9) =sin 2 (6»)/4. 
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4 Fourth-order Gram-Charlier expansion series 

In case of a detailed transition array, the Zeeman broadening of a line can be represented by a fourth-order 
A-type Gram-Charlier expansion series: 



~i » B By/2*MM V 2 V 3 / 
+^^(3-6 y? 2 + 4)), 



where y g = (E — E^j n >ji — q g e (abB) / '((ibB-J M-^i the coefficient <7 e of the line 7J — > 7' J' being 
defined in equation (U]). If the values of <? 7 j and <jyj' are unknown, we suggest to replace g e by its average 
value in LS coupling <? e , which can be roughly estimated, knowing the distribution of spectroscopic terms 
Q(S,L) [Gilleron & Pain (2009)], by: 

9e=J2 E Q(S,L')g e (S,L,J,L',J') A(L,L',J,J'), (9) 

S,L,J L',J< 
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Figure 1: Effect of a 1 MG magnetic field on triplet transition ls2s 3 S — > ls2p 3 P of C V. The Doppler 
full width at half maximum (FWHM) is taken to be 0.005 eV. Since, for a Gaussian, FWHM=2 yjl ln(2)u, 
this corresponds to a variance v ~ 4.5 10 -6 eV 2 . The observation angle 8 is such that cos 2 (6') = 1/3. 



where A(L, L' , J, J') stands for the selection rules: \L — 1| < L' < L + 1 avoiding V = L = and 
| J — 1\ < J' < J + 1 avoiding J' = J = 0. The notation g e (S, L, J, L', J') means g e in which the Lande 
factors are estimated in LS coupling: 

3 S(S + 1)-L(L+1) 
5S "=2 + 2J(J+1) ' (10) 

Figure [T] shows, in the case of carbon 4+ (C V in spectroscopic notation) for a 1 MG magnetic field, 
that Gram-Charlier series ((8| provides a good depiction (see the black curve named "Statistical" in the 
legend) of the triplet transition ls2s 3 S — > ls2p 3 P. 

5 Conclusion 

We proposed a statistical modelling of Zeeman patterns that can be easily implemented in opacity codes, 
together with a simple procedure for the diagnostic of magnetic fields. 
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